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Abstract 

Electric signals have been recently recorded at the Earth's surface with amplitudes appreciably 
larger than those hitherto reported. Their entropy in natural time is smaller than that, S u , of 
a "uniform" distribution. The same holds for their entropy upon time-reversal. This behavior, 
as supported by numerical simulations in fBm time series and in an on-off intermittency model, 
stems from infinitely ranged long range temporal correlations and hence these signals are probably 
Seismic Electric Signals (critical dynamics). The entropy fluctuations are found to increase upon 
approaching bursting, which reminds the behavior identifying sudden cardiac death individuals 
when analysing their electrocardiograms. 
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The time series analysis of various phenomena in complex systems (and especially those 
associated with impending catastrophic events, e.g. Q, Q|) in the framework of the newly 
defined time-domain [ 
ectrocardiograms P, 



, ^], termed natural time, reveals interesting features. Examples are 
seismicity[3] and Seismic Electric Sig nals (SES) activities^, U, fl 
Si y, llQ| • This new time domain is optimal for enhancing the signals in time-frequency space 
when employing the Wigner function and measuring its localization property 11 1. 

In a time series comprising iV pulses, the natural time Xk = k/N serves as the index ^, 0] 
for the occurrence of the k-th event. It is, therefore, smaller than, or equal to, unity. 
In natural time analysis, the time evolution of the pair of the two quantities (xk, Qk) is 
considered, where Qk denotes the duration of the fc-th pulse. The entropy S in the natural 
time-domain is defined 2| as the derivative of the function (x q ) — (x) q i n respect to q, for 



q = 1, which gives |3,la|: 



S= (xlnx) - (x)ln(x) 



(1) 



where (f(x)) = J2k=iPkf(Xk) and 
on the sequential order of pulsesp, ll 



? fc = Qk/ ^2n=fi,Qn- It is a dynamic entropy depending 
] and exhibits^] concavity, positivity and LescheQ, Q 



stability. When the system enters the critical stage (infinitely ranged long range temporal 



correlations the S-value is smaller^, |10| than the value $„(■■ I ,/2 hi 2 - I / I O.OOGO) 
of a "uniform" distribution (defined in Refs. P, U 
i.e., 

S < S u 



e.g. when when all pk are equal), 

(2) 



The value of the entropy obtained upon considering the time reversal T, i.e., Tp k =_PN-k+i, 
is labelled by S-. An important point emerged from the data analysis in Ref.^l is the 
following: Although the study of the S values enables the distinction between SES activities 
and noises produced by nearby operating artificial (man made) electromagnetic sources 
(AN), i.e., S < S u for the SES activities and S ^ S u for AN, this does not hold for the S_ 
values. This is so, because for the SES activities we found that the S- values are smaller 
than (or equal to) S u , while for AN the S- values are either smaller or larger than S u . Here, 
we provide more recent data on the SES activities, which strengthen the conclusion that 
both S and S_ are smaller than S u . In other words, the following key point seems to hold: In 
signals that exhibit critical dynamics (i.e., SES activities) upon time- reversal their entropy 
values (though become different than those in forward time) still remain smaller than S u . 
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Why? The answer to this question is a challenge, because, if it is generally so, among 
similar looking signals we can distinguish those that exhibit critical dynamics. Since the 
latter signals are expected to exhibit infinitely ranged long range correlations, this might be 
the origin of the aforementioned behavior. To investigate this point, numerical simulations 
are presented here for fractional Brownian motion (fBm) time series as well as for an on-off 
intermittency model (for the selection of the latter see Ref . ) . The simple case of an fBm 
was selected in view of the following suggestion forwarded in Ref.j?| concerning the Hurst 
exponent H . If we assume that, in general, H is actually a measure of the intensity of long 
range dependence, we may understand why in SES activities, when analyzed in the natural 
time domain, we find H values close to unity, while in AN (where the long range correlations 
are weaker^) the if-values are markedly smaller (e.g., around 0.7 to 0.8). 

We first present the recent experimental results. Figure 1 depicts five electric signals, 
labelled I>. M 3 , M 4 and V\, that have been recorded on March 21, 23 and April 7, 
2005. Note that the first four (Mi to M4) have amplitudes that not only are one order 
of magnitude lar ger than the fifth one (Vi, see Fig.lc), but also significantly exceed those 
hitherto reported (l5j. The way they are read in natural time can be seen in Fig. 2. The 
analysis of these signals leads to the S and 5*_ values given in Table HI an inspection of which 
reveals that they are actually smaller than S u . Hence, on the basis of the aforementioned 
criterion, these signals cannot be classified as AN, but they are likely to be SES activities. 
Note that, although in general S is different than 5_, no definite conclusion can be drawn 
on the sign of S — S- (see also Table I of Ref.p). 

We now present our results on fBm. We first clarify that Weron et al. recently studied 
an algorithm distinguishing between the origins (i.e., the memory and the tails of the process) 
of the self-similarity of a given time series on the base of the computer test suggested in 



Ref. [Uj. By applying it to the SES activities, they found the fBm as the appropriate type 
of modelling process. The fBm, which is if -self similar with stationary increments and it is 
the only Gaussian process with such properties for < H < can be simulated jlfl I20I] . 

see also pp.321 -323 of Ref.Q, by randomizing a construction due to Weierstrass , i.e., using 



the Weierstrass-Mandelbrot function 
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^ sin(fr ; t + di) 
1 

where b > 1, q normally distributed with mean and standard deviation 1, and d\ are 



uniformly distributed in the interval [0, 2tt] (cf. when using the increments of Eq.fjSJ) one 
can also produce fractional Gaussian noise of a given H). By using Eq. (j3J), ffim for various 
values of H were produced, the one-signed segments of which were analyzed in the natural 
time domain (an example is given in Ref.Q]). The Monte-Carlo simulation results for each 
segment include not only the values of the entropies S and S"_, but also the exponent «dfa 
of the Detrended Fluctuation Analysis (DFA)[3, [^J. For segments of small number of 
points iV (cf. only segments with iV > 40 were considered), the value of oidfa may vary 
significanly from that expected for a given value of H; DFA was preferred, because it is one 
of the few well-defined and robust estimators of the scaling properties for such segments(e.g. 
0], see also pp. 300-301 of Ref.jnJ). The results are shown in Fig.3, in which we plot 
the S and S_ values versus Q-dfa- Since the analysis of the SES activities in natural time 
result 0, Q in DFA exponents cxdfa around unity, hereafter we are solely focused in Fig.3 
in the range 0.8 < «dfa < 1-2. An inspection of this figure reveals the following three 
conclusions: First, despite the large standard deviation, we may say that both S and S- are 
smaller than S u (~ 0.0966) when cxdfa ~ 1- Second, S and S- are more or less comparable. 
Third, comparing the computed S and S- values (~ 0.08 for aoFA ~ 1) with those resulting 
from the analysis of the SES activities (see Table HI see also Table I of Ref. 9]), we find a 
reasonable agreement. Note that these computations do not result in a definite sign for 
S — S_ in a similar fashion with the experimental results. 

In what remains, we present our results on a simple on-off intermittency model. We clarify 
that on-off intermittency is a phase-space mechanism that allows dynamical systems to 
undergo bursting (bursting is a phenomenon in which episodes of high activity are alternated 
with periods of inactivity likewise in Fig.2e). This mechanism is different from the well 
known Pomeau-Manneville scenario for the behavior of a system in the proximity of a saddle- 

n 

node bifurcation [25]. Here, we use the simple model of the driven logistic map 

X t+1 = A(Y t )X t [l - X t ] (4) 

where we assume that the quantity A(Y t ) is monotonic function of Y t and that < A < 4 
(cf. A is further specified below). The system has the invariant manifold X = and the 
level of its activity is measured by Atj^]. In order to have the on-off mechanism in action, 
we specialize to the case of a noise-driven logistic map, with 

A(Y t ) =A + aY t (5) 



where Y t is (^-correlated noise which is uniformly-distributed in the interval [0,1] and Aq and 
a are parameters. In order to have < A < 4, we assumej^G] Aq > 0, a > and Aq + a < 4. 
The relevant parameter plane for the noise-driven system of Eqs.(@J) and (j5J) (as well as the 



parameter range for which the fixed point X = is stable) can be found in 



28, 



ig. 1 of Ref 



291 ]. Bursting is 



while the description of the intermittent dynamics is given in Refs.[27, 
observed in the temporal evolution of X t as the stability of the fixed point X = varies. 
Following Ref.j3], for Aq = there is a critical value a c > 1, below which the system 
asymptotically tends to the fixed point X = 0, without any sustained intermittent burs ting . 
For this case, i.e., Aq = 0, the value a c = e = 2.71828. . . leads to on-off intermittency 26] . 
In the intermittent system under discusssion, both the signal amplitude and the power 
spectrum resulted j3] in power-law distributions (with low frequencies predominating in the 
power-spectrum) 

Several time-series have been produced for the above on-off intermittency model with 
the following procedure: The system was initiated at a time (ti n = —200) with a uniformly 
distributed value X tin in the region [0,1], and then the mapping of Eqs.(0J) and (jnj) was 
followed until N events will occur after t = 0. The results for X t ,t = 1,2... AT were 
analyzed in natural time domain and the values of S and S_ have been determined. This 
was repeated 10 3 times (for a given number, N, of events) and the average values of S and 

have been deduced. These values are plotted in Fig.4(a) versus (a — e)N 1 / 2 (The factor 
iV 1 / 2 stems from finite size scaling effects, since for large values of N, e.g., iV >15000, a 
scaling -reminiscent of a 1st order phase transition- was observed, details on which will be 
published elsewhere). This figure reveals that as the critical value is approached from below 
(i.e., a — ► e_) both S and are smaller than S u . Note that Fig.4(a) also indicates that S 
is probably larger than 5 1 -, while in the fBm time series no definite sign for S — S- could 
be obtained. 

Another interesting point resulted from the on-off intermittency model is depicted in 
Fig. 4(b), where we plot the fluctuations 5S and 5S- (i.e., the standard deviations of the 
entropies S and S-, respectively) versus (a — e)N 1 / 2 . It is clearly seen that these fluctu- 
ations are dramatically enhanced as the critical value is approached (i.e., a. — > e). This is 
strikingly reminiscent of our earlier results 0, |(| upon analyzing electrocardiograms (ECG) 
in natural time domain and studying the so called QT intervals. These results showed that 
the fluctuations of the entropy SS(QT) are appreciably larger in sudden cardiac death (SD) 
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individuals than those in truly healthy (H) humans (see Fig. 2 of Ref.[6j). We emphasize 
that the aforementioned points should not be misinterpreted as stating that the simple logis- 
tic map model treated here can capture the complex heart dynamics, but only can be seen 
in the following frame: Since sudden cardiac arrest (which may occur even if the electrocar- 
diogram looks similar to that of H) may be considered as a dynamic phase transition^y], it 
is reasonable to expect that the entropy fluctuations significantly increase upon approaching 
the transition. 

In summary, recently recorded electric signals (having the largest amplitudes recorded 
to date) exhibit the property that both S and S- are smaller than S u and hence are likely 
to be SES activities (critical dynamics). This property seems to stem from their infinitely 
ranged long range correlations as supported by computational results in: (1) fBm time 
series and (2) a simple on-off intermittency model. The latter model also suggests that the 
fluctuations (SS and 5S_) significantly increase upon approaching the transition, which is 
strikingly reminiscent of the increased (^-values found for the QT-intervals for the sudden 
cardiac death individuals. 
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FIG. 1: Electric Signals recorded on March 21, 2005 (a), March 23, 2005 (b) and April 7, 2005 (c). 
The signals in (a) and (c) are labeled hereafter Mi and V\ respectively, while that in (b) consists 
of the three signals' activities labeled M2, M3, M4 (sampling frequency / exp =lHz). The Universal 
Time (UT) is marked on the horizontal axis. Additional details for the two dipoles -records of 
which are shown here- as well as for the sites of the measuring stations are provided in Ref. \l\ . 
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FIG. 2: How the signals depicted in Fig.l are read in natural time, (a), (b), (c), (d), (e), correspond 
to the signals' activities labeled Mi, M2, M3, M4 and V\, respectively. 



TABLE I: The values of S and S- together with the number of pulses iV for the SES activities (the 
original time series have lengths between 2 x 10 3 and 10 4 , compare Fig.l with f exp =lRz) shown in 



Fig.l. 

Signal N S S- 

Mi 78 ±9 0.094±0.005 0.078±0.003 

M 2 103 ±5 0.089±0.003 0.084±0.003 

M 3 53 ±3 0.089±0.004 0.093±0.004 

M 4 95 ±3 0.080±0.005 0.086±0.006 

Vi 119 ± 14 0.078±0.006 0.092±0.005 
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FIG. 3: Calculated values of S (squares) and 5_ (triangles) versus the DFA exponent odfa- The 
error bars indicate the region of one standard deviation (±cr). The horizontal line corresponds to 
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FIG. 4: (Color Online) Calculated results for the on-off intermittency model discussed in the text: 
(a) depicts the average values of S (closed symbols) and S- (open symbols) while (b) those of the 
fluctuations 5S and 6S- versus the finite size scaling variable (a — a^N 1 / 2 . The quantity iV stands 
for the number of the events considered in each sample time series; iV=70000, 50000, 30000, 15000 
correspond to squares, circles, triangles and inverted triangles, respectively. The horizontal line in 
(a) corresponds to S u . 
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